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Circulations 
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Transhorizon vhf and uhf fields exhibit deep fades or large signal enhancements of several 
hours' duration, as the propagation mechanism alternates between partial reflection and 
scattering caused by turbulent dielectric fluctuations in the atmosphere. Such alternations 
occur when strong refractive layers develop below 3,000 ft. Surface wind streamline analyses 
show that mesoscale centers of convergence or divergence cause local redistribution of 
refractive layering, tending to produce the change from one mechanism to the other. 

Current scattering theory and the empirical findings of others are examined to determine 
the gross meteorological factors that influence changes in scattered fields. The two variables 
in the turbulent scattering coefficients, the scattering angle and the intensity of dielectric 
fluctuations at high wave numbers, are found to be dependent upon the refractive layering 
and the thermal stability of the airmass. It has been shown elsewhere that refractivity 
and stability are principally functions of the vertical velocity in the atmosphere. It is shown 
here that the direction and relative magnitude of the vertical velocity can be inferred from 
the upper-tropospheric wind velocity divergence. Received scattered signals are found to 
be well correlated with computed velocity divergence. 

It is suggested that the variations of scattered signal level or range can be predicted in 
a routine manner by regular meteorological personnel using ordinarily available meteoro- 
logical data. 



1. Introduction 

Future Navy communication networks and early 
warning radar may depend to a great extent upon the 
tropospheric scatter mode of radio propagation. 
The development of high-gain low-noise receivers 
and statistical target identification techniques makes 
the vhf and uhf region of the radio spectrum useful 
for ranges well beyond the radio horizon. 

The signal in the far transhorizon region is char- 
acterized by a rapidly fluctuating power amplitude 
superimposed on a more slowly changing average 
signal level. The rate of fast fading depends on 
distance and radiofrequency, indicating the existence 
of {uii correlated multipath transmission from num- 
erous moving scattering centers. The slow fading 
has seasonal, average diurnal, and random day-to- 
day or hourly components. Seasonal and average 
diurnal fluctuations are amenable to statistical 
analysis and are predictable within relatively nar- 
row limits) [l]. 2 The random daily or hourly varia- 
tions in signal level have been less well understood. 
These are the fluctuations in signal level (or range) 
in which the Navy may be most interested. 

The best current theories explain the scatter propa- 
gation mechanism by making use of scattering co- 
efficients which are functionally related to the gross 
refractivity of the lower atmosphere and the intensity 
of turbulent dielectric fluctuations at radiofrequency- 
dependent scale sizes. The scattering coefficients 
are difficult to measure and, at the present time, are 
unpredictable. It is possible to predict some aspects 
of the atmosphere's refractive structure which, from 
elementary considerations of theory, are most likelv 
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to control the magnitude of the scattering coeffi- 
cients and hence the transhorizon signal strength. 

This report indicates the meteorological factors 
which appear to control the electromagnetic scatter- 
ing mechanism and suggests an operational method 
for predicting the direction and relative magnitude 
of signal level change. 

2. Discussion of Theoretical and Empirical 
Relationships 

The theories proposed to explain the large other- 
than-diffraction transhorizon fields are based on two 
propagation mechanisms: (1) partial reflection of in- 
cident radiation by horizontally stratified dielectric 
layers [2, 3, 4], and (2) scattering by turbulent dis- 
continuities in the dielectric structure of the atmos- 
phere [5, 6, 7]. The phenomenon of partial reflection 
has been recognized and studied for over ten years. 
In the microwave region, for usual refractive gradi- 
ents, the effect of layers decreases exponentially with 
height, so that reflections of this sort become unim- 
portant for layers above about 3,000 ft. More 
recently, Bauer, Friis, et al. [8, 9] theorized that 
energy may be partially reflected in much the same 
manner from thin, highly refractive sublamina within 
the region illuminated by both antennas. Refrac- 
tometer measurements made during convective at- 
mospheric conditions by Navy Electronics Labora- 
tory (NEL) failed to show the existence of such 
layering at times when transhorizon radio transmis- 
sion was normal. It is likely that, under conditions 
of a well-established elevated layer, both propagation 
mechanisms may exist concurrently [10], but that in 
the absence of strong relatively low elevated layers, 
the turbulent scattering mechanism is dominant. 
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The turbulent scattering theories, following the 
original work of Booker and Gordon [5], postulate in 
common an equation of the general form: 



Proz- 



F(K). 



(1) 



sin" 



Pr= received power, 
0= intersection angle of rays from transmitter 
and receiver, 

K=4tTT sin - 



X= radio wavelength, 
F(K)= function of dielectric fluctuation intensity at 
wave number, k, and 
n= exponent, dependent on form of F. 

By incorporating the modern theories of turbulence 
put forth by Obukhoff, Kolmogoroff, and others with 
radio scattering theory, Batchelor [7] found that, for 
an isotropic radiator, n has the value 11/3. For 
propagation with narrow beam antennas, the expo- 
nent becomes 14/3. It should be noted that for 
grazing propagation and ducting gradients, expres- 
sion (1) fails; propagation is line-of-sight and, dis- 
regarding focusing effects, power is maximized. 

2.1. Scattering Angle 

The strong dependence of signal level upon scat- 
tering angle (mean refractive gradients) is demon- 
strated in figure 1 . The abscissa is the signal level 
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Figure 1. Scattered field signal level as a function of scattering 
angle and refractive gradients in a homogeneous refractive 
atmosphere normalized to standard gradient. 



in decibels above that received in a standard atmos- 
phere (— 12 N units/ 10 3 ft) for ratios of the standard 
atmosphere scattering angle to the actual scattering 
angle (ordinate on the left) given by the N gradients 
(along the right of the graph). It is noteworthy 
that with high signal levels and near-ducting gra- 
dients, deep fades usually associated with other 
modes of propagation can occur with very slight 
changes in the low level refractive gradients. 

A method for determining the scattering angle 
from the refractive index profiles over the scatter 
propagation link is outlined in appendix A. The 
method is generalized for a multilayered atmosphere, 
the layers having a limited degree of horizontal 
homogeneity. 

Other investigators concerned with developing 
methods for predicting transhorizon fields through 
use of meterorological parameters have found cor- 
relation between scattered fields and low-level re- 
fractive gradients. Bean [1], reasoning that lower 
tropospheric refractive gradients are directly related 
to the surface value of refractive index, found 
statistical correlation between transhorizon fields 
and the surface refractivity. Misme [11], using a 
weighted average tropospheric refractivity, found 
correlation between this "useful gradient" and the 
received signal. 

Time-averaged low-level refractive gradients have 
been shown by Moler and Arvola [12], to be princi- 
pally a function of the large-scale vertical motions 
in the atmosphere. Later in this report, computa- 
tions of convergence at the 300-millibar (mb) level 
and the implied subsidence (see app. B) throughout 
the troposphere are compared with the development 
and descent of a strong refractive layer and the 
resultant change in scattering angle and radio signal 
strength. 

Low-level gradients may also be temporarily 
altered with the development of minor eddies in 
the surface wind flow. Occasionally, during the 
NEL scatter propagation experiments, radio signal 
levels and characteristics would indicate a change 
in the propagation mechanism from scattering to 
partial reflections. Subsequent analysis of the small- 
scale wind circulations showed that local regions of 
convergence (divergence) which develop within a 
relatively homogeneous large-scale weather regime 
cause raising (lowering) of refractive layers and, 
consequently, a change from one dominant mode of 
propagation to another. Examples of this will be 
given in a later section of this report. 

2.2 Turbulent Dielectric Fluctuations 

Little is known about the variation of intensity 
of turbulent dielectric fluctuations at wave numbers 
important for microwave scattering as a function 
of the meterorological situation. Recent work by 
Gossard [13] indicated that the refractive index 
spectral intensity at high wavenumbers in the inertial 
subrange of turbulence is only weakly dependent 
upon atmospheric stability or refractive gradients. 
From this work, it may be concluded that variations 
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in signal level result almost entirely from refractive 
changes in scattering angle or that there are thin 
zones of quite variable turbulent intensity which 
superimpose an amplitude modulation upon a low 
average signal. 

Batchelor [7] suggested that the fluctuations of 
refractive index may be due to turbulent convection 
in the presence of a spatial gradient of mean refrac- 
tive index. If this suggestion is correct, it seems 
likely that the greatest intensity of dielectric fluctua- 
tion would occur along the lower boundary of 
elevated stable layers (refractive layers) where con- 
vective energy from the lower, less stable region is 
absorbed and dissipated. 

Based on empirical evidence, Misme [11] found 
that "for equal refractive index gradients, the atmos- 
phere which is characterized by the greatest stability 
permits the highest fields to be received from a 
great distance. " By "equal gradients", he meant the 
refractivity which determines the scattering angle. 

Gordon [14], using assumed values of fluctuation 
intensity, concluded that the contribution of scat- 
tered energy from a single turbulent layer may be 
greater than that from the remaining common 
volume. The theories of partial reflection [2, 3, 4, 
8, 9] require the presence of highly refractive sub- 
strata within a mean refractive gradient with greatest 
radio fields associated with intense elevated refrac- 
tive layers. 

All evidence indicates that departures of radio 
fields from a low ambient level can best be forecast 
indirectly through a prediction of the atmospheric 
thermal stability. It is well known that, in general, 
atmospheric stability is largely a function of the 
vertical components of the winds. At the onset of 
subsidence, refractive layers develop; continuous 
subsidence causes lowering and intensification of the 
layers and increases the thermal stability of the air- 
mass. Rising motion tends first to elevate the layers 
and later to cause their complete dissipation. 

The behavior of radio signal levels under the in- 
fluence of positive as opposed to negative vertical 
motion will be shown in a later section. The method 
for deducing the sign and relative magnitude of the 
vertical wind component from upper tropospheric 
constant pressure maps will be found in appendix B. 

3. Operations 

In January 1958, NEL activated three trans- 
horizon radio links which were in periodic operation 
until mid-April. Stations at Santa Barbara (190 
statute miles) and at Point Mugu (144 statute miles) 
transmitted 9,380 Mc (X-band) and 1,365 Mc (Z- 
band) pulsed signals. A third station on San 
Clemente Island (78 statute miles) transmitted at 
9,380 Mc only. All receivers were located on Point 
Loma, San Diego. Figure 2 is a map of the southern 
California operating area. Signal levels used in this 
report were measured over 1-min intervals on 
integrating equipment and recorded in decibels with 
respect to 1 milliwatt. 
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Figure 2. Navy Electronics Laboratory scatter propagation 
experiment operations area. 



Strong low-level elevated refractive layers are a 
dominant feature along the southern California coast 
and frequently cause ducting or abnormally strong 
signal enhancement from internal reflections when 
the layers are lower than about 3,000 ft [3, 4]. In 
order to avoid such conditions, the links were placed 
in operation only when there was a reaonable assur- 
ance that refractive layers would be absent or at a 
high elevation. Synoptic experience has shown that 
airmasses are most unstable and thermal layering 
absent when under the influence of a strong cyclonic 
circulation through most of the troposphere. For 
this reason, the bulk of data was taken during periods 
when a deep low pressure system at the 500-mb level 
(18,000 ft) was near or over the radio link. 

Radio link operations were discontinued whenever 
signal characteristics indicated that a propagation 
mode other than scattering was dominant and likely 
to continue. Unfortunately, as a result very little 
data were taken during periods of high signal level 
and conditions of airmass subsidence. 

The southern California sectional maps were con- 
structed from hourly airways teletype reports and 
from weather logs of two aircraft carriers which were 
operating in coastal waters. The surface synoptic 
and 500-mb charts are copies of National Weather 
Analysis Center maps transmitted over the facsimile 
network. Calculations of vorticity advection at the 
300-mb levels were made from maps plotted at NEL, 
using data from the Daily Bulletin of Northern 
Hemisphere Data Tabulations published by the 
U.S. Weather Bureau. 

4. Measurements 

4.1. Signal Level and Refractive Layering 

The constant pressure contours for the 300-mb 
pressure surface for 1600 P.s.t., March 16 to 19, are 
shown in figure 3. The field of relative vorticity 
(dashed lines) indicates positive vorticity advection 
(horizontal divergence) east of the upper level trough 
and negative vorticity advection (horizontal con- 
vergence) west of the trough. From continuity con- 
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Figure 3. Three hundred millibar pressure-height contours (solid) and isopleths of relative vorticity (dashed) in units of 10 5 t 
at 1600 P.s.t. for (a) 16 March (b) 17 March (c) 18 March (d) 19 March. 



siderations (app. B), positive vorticity advection 
implies negative vertical velocities. Figure 4 shows 
V'V 300 mb as a function of time compared with a 
series of refractive index profiles derived from the 
Los Angeles radiosonde observations. The func- 
tional relationship between the scattering angle 6 
and the refractive gradients has been derived in 
appendix A. Figure 5 shows the scattering angle for 
the grazing ray and height of the common volume for 
the Santa Barbara link during the 3-day period. 
The upper curve shows the signal level as a function 
of the scattering angle normalized for a standard 



atmosphere. The second curve is the radio field 
strength for the Santa Barbara A'- band link for the 
period 0800 to 1600 P.s.t,, 19 March. Radio data 
were not taken the previous days. 

On March 26 there was an anomalous 20-db 
enhancement of signal level during a 3-hr period near 
noon that appeared unexplainable on the basis of 
radiosonde data or by vertical velocities inferred 
from the upper-troposphere divergence field. The 
refractive index profiles reduced from Los Angeles 
radiosonde observations for 0400 P.s.t, and 1600 
P.s.t. show an elevated layer near 4,000 ft (fig. 6) 
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Ficure 4. Hoi izontal velocity divergence at the 300-mb pressure 
level (upper graph) and Los Angeles and San Diego refractive 
index profiles for the period 16 to 1!) March 1958. 
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Figure 5. Upper curve indicates power depending upon 
scattering angle, 6, for the period 17 to 19 March 1958. 
Second curve shows the received power in decibels per milliwatt 
for Santa Barbara X band. The lower two curves show the 
height of the scattering volu?ne, h v , and the scattering angle, 6, 
during the same period. 

and computations of V-Fmade at 300 mb (upper 
curve fig. 6) indicate a continuous reduction in 
stability taking place. A series of isobaric and 
streamline southern California sectional charts were 
analyzed at hourly intervals for the day. Figure 7 
shows the map series and includes the large-scale I 
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FIGURE (>. Refractive index profiles nil need from Los Angele 
rabiosonde <in<l (upper curve) 300-mb velocity divergence for 2( 
to 27 March 1958. 



synoptic chart for 1000 P.s.t. 



8 is the Santa 



From figure 7, 



igun 
Barbara X-band signal lor the day. 
at 0900 P.s.t., a cyclonic eddy with two centers 
northwest of San Diego lav across the Santa Barbara 
radio link. Streamlines indicate convergence and 
consequent upwelling of the surface air in the con- 
vergent areas. By 1 100 P.s.t. the cyclonic (h\(\\ had 
weakened, due to \\\r onset of a sea breeze and" sub- 
sequent low-level divergence and subsidence of air 
over the Catalina channel. The 1200 P.s.t, map 
shows only a weak remnant of a zone of convergence 
south of Point Mugu and signals have risen 20 db. 
By 1300 P.s.t. offshore convergence had ceased 
entirely and signals had reached a peak level for the 
day shortly before 1300 P.s.t. The 1300 P.s.t. map 
indicates a strengthening of the pressure trough north 
of Point Mugu with a tendency for the trough to 
extend seaward. Signal levels started to fall after 
1300 P.s.t. By 1500 P.s.t., the convergence zone 
had undergone strong redevelopment across the 
Santa Barbara link and signal levels had fallen about 
22 db from the peak level. 

4.2. Signal Level and Upper-Tropospheric 
Divergence 

During the period 26-28 March inclusive, the 
radio links were in operation from 0800 P.s.t. to 1600 
P.s.t. daily. The weather regime was dominated by 
an upper-level cyclone which developed off the 
northwest coast of California and moved southeast- 
ward across the radio link. A series of 500-mb 
height contour charts drawn at 12-hr intervals during 
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Figuke 8. Santa Barbara X.-band signal level for 26 March 

1958. 



the period are shown in figure 9. The isopleths are 
contours of the height of the 500-mb pressure sur- 
face drawn for 200-ft intervals. The velocity 
divergence along the streamlines 5° upstream from 
the operating area was computed for each map dur- 
ing the period. Figure 10 is a graph of the velocity 
divergence at 500 mb and average signal level re- 
ceived from the Santa Barbara i-band and Point 
Mugu X-band transmitters, both plotted as a func- 
tion of time. The signal levels are average values 
for the periods 0800-1000 P.s.t., and 1400 to 1000 
P.s.t. plotted at 0900 and 1500 P.s.t., respectively, 
and connected with a smooth curve. Solid portions 
ofjjthe curve are for hours of radio transmission. The 
dashed curves are extrapolations of the solid curve 
slopes. Unfortunately, periods of radio silence were 
longer than transmission periods so that one can 
only assume that the signal strength curves are 
reasonably close to those extrapolated. Although it 
would be impossible to predict absolute signal level, 
the correlation between the sign and relative magni- 
tudes of the time derivatives of the curves is good. 

Computations of vorticity advection were made 
at the 500-mb pressure level for all days when there 
was a reasonably good sample of radio data and at 
the 300-mb level for periods of particular meteoro- 
logical interest. Signal levels for the 2-hr period 
nearest map time were averaged for the conditions of 
positive and negative vorticity advection (positive 
and negative vertical velocities). Table 1 gives these 



Table 1. Average values of signal level for conditions of positive 
and negative vertical wind velocities 
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-67.0 
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-57.6 
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-62.3 
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-60. 


M x 
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Ml.. -- 


-64.4 







average values of signal under the two conditions 
from a sample size of _33 at the 500-mb level and 18 
at the 300-mb level 



Both the X- and L-band signal 



levels are significantly higher under conditions of 
subsidence with no apparent frequency dependency 
to be noted. 

5. Results and Conclusions 

Other surface streamline charts and upper-tropo- 
spheric synoptic series were analyzed, but to presenl 
the results here would be repetitious and redundant. 
All concurrent meteorological and radio signal level 
analyses showed essentially the same features as 
indicated in the previous sections. 

During the NEL experiment, sea level values of 
refractive index were essentially invariant with time, 
yet signal levels ranged over approximately 60 db. 
It seems clear that for over-water paths, short-term 
propagation characteristics are independent of sur- 
face refractivity. 

Nonducting low-level refractive layers influence 
transhorizon fields in two ways: (1) refraction and 
subsequent decrease of scattering angle, and (2) 
partial reflection at the interfaces of superrefractive 
layers as the layers approach critical reflection 
heights and intensities. The second effect dominates 
for layers below about 2,000 ft, For example, if after 
Batchelor [7], a dielectric fluctuation spectrum of the 
form E (K) oc R~$ is accepted, then for narrow-beam 
propagation, the received power dependency upon 
scattering angle would be Pr oc 0/2~~s\ Assuming 
standard propagation early in the morning of March 
19, the Los Angeles sounding for 1600 P.s.t. March 19 
(fig. 4) would indicate a lessening of scattering angle 
sufficient to cause a 12.7-db rise in signal level (fig. 
5). The actual rise in signal was approximately 35 
db. In order to account for such a signal rise by 
refraction, the scattering angle would have had to 
decrease by nearly one-sixth. It is unlikely that the 
overwater refractive gradients were sufficiently 
greater than those shown by the Los Angeles 
radiosonde to account for this much refraction. 
One is forced to conclude that the partial reflection 
mechanism becomes more important for transhorizon 
propagation than scattering when relatively low 
layers are present. This suggests that the refractive 
index climatology charts for oceanic areas [15] 
could be useful for determining the distribution of 
dominant modes of transhorizon propagation. 

The short term prediction of change of radio range 
for oceanic areas becomes greatly simplified if the 
problem can be looked upon as one concerning only 
changes in thermal stability. Except for coastal re- 
gions where continental air may undergo rapid mod- 
ification, refractive layering and thermal stability 
over the oceans are principally functions of vertical 
wind velocities. Techniques for computing exactly 
or making qualitative estimates of the vertical wind 
velocities similar to the one outlined in appendix B 
are familiar to meteorologists. It is likely that 
satisfactory estimates of radio scatter propagation 
ranges can be made indirectly from regularly avail- 
able weather charts. 
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Figure 9. Five hundred millibar 'pressure-surface contour charts at 12-hr intervals from 0400 P.s.t. 26 March to 0400 P.s.t. 

28 March. 
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Figure 9. — Continued. 
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Figure 10. Santa Barbara L-band and Point Mugu X-band 
signal levels and 500-mb velocity divergence for 0400 P.s.t. 
for\26 March to 1600 P.s.t. 28 March inclusive. 
Dashed signal level curves are extrapolated values for periods of radio silence. 



6. Appendix A: Scattering Angle and Height 
of the Common Scattering Volume 

The scattering volume is the portion of atmosphere 
common to both transmitting and receiving antenna 
beams. The acute angle at the beam intersection is 
called the scattering angle. It is possible to com- 
pute the height of the base of the scattering volume 
and the magnitude of the scattering angle from the 
gross refractive index profile reduced from radiosonde 
or refractometer data. 



The atmospheric model (fig. 11) assumes horizon- 
tal homogeneity of a refractive layer for the distance 
that a ray segment remains within the layer. The 
number of refractive layers may vary with distance 
over the scattering link. The horizontal distribution 
of vertical refractive index gradients may be inter- 
polated from soundings taken near either end of the 
radio link. 

In figure 11, unprimed quantities refer to param- 
eters between the transmitting (receiving) antenna 
and the beam intersection; primed quantities refer 
to parameters between the intersection and the 
receiving (transmitting) antenna. Symbol meanings 
are as follows: 

A= antenna heights, 

<i = distance to the radio horizon, 

d t = projection on the earth of the ray path in 

the i-th layer, 
D= distance between antennas, 
Wi= thickness of i-th layer, 
M= modified refractive index = (n— 1) X 10 6 

-f 0.048 h, 
h v = height of base of the common volume, 
0= scattering angle. 

Subscript v refers to the refractive layer wherein 
the beams intersect. From figure 11: 



S=d v +d' v =D-Jl d { . 



(Al) 



The angle which the ray makes with the local 
horizontal upon emerging from the n-th layer is given 
by 



[■ 



7 { ^iJ. 



1 



dh 



(A2) 



For the case of grazing propagation a is the angle 
the ray makes with the horizon and is zero. For 
nongrazing propagation a Q is the elevation angle of 
the antenna beam. Figure 12 is a graph of 2WdM/dh 
as a function of w and AM/1,000 ft. Figure 13 is a 
nomographic solution of (A2). The incremental 
distances d Q and d^ can be easily determined from 
figure 14, the solution of 

fci °'Hfy'T- <A3 > 

The parameter H is either antenna height or the 
depth of the surface-based refractive layer at each 
end of the link. 

The incremental distances d t from d 2 to d v -\ may be 
determined from figure 15, a graphical solution of 



d i =10 & (a i —a i 



-Km 



(A4) 



S may now be determined according to (Al). 
The remaining distances, d„ and d,' , are given by 
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and 



* ?(f).'°-+^-+^ 

1(f). lo-^-W-, 



d v =S — (i . 



(A5) 



(A6) 



The height of the base of the common volume is now 
determined from 

and the scattering angle from 



ah „ 



(A8) 




Figure 11. Schematic diagram of scatter propagation link with nonhomogeneous, layered 

atmosphere. 
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Figure 13. Nomogram of a n =f(Wi, AM/10 3 ft, a ). 
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Figure 14. Nomogram of d = f(H, AM/lCPit). 



7. Appendix B 

7.1. Method for Estimating the Relative Magnitude 
and Direction of the Vertical Component of the 
Wind 

It may become necessary for the radio meteorolo- 
gist to prepare operational propagation forecasts for 
areas where machine-computed vertical velocity 
charts are not available. The following simple 
method for obtaining a "better-than-guess" estimate 
of the vertical -mot ion field is suggested as a forecast 
aid . 

The atmospheric model assumes the following 
conditions: 

1 . An equivalent barotropic level (or level of non- 
divergence) exists in the mid troposphere, usually at 
a lower altitude than the 500-mb surface and cer- 
tainly below the 300-mb surface. 

2. The horizontal wind velocity divergence changes 
sign at the level of nondivergence (from here on 
"LND") on the vertical coordinate and at no other 
elevation. 

3. The direction and magnitude of vertical velocity 
throughout the troposphere are the same in sign and 
proportional in magnitude to those at the LND. 

4. Vertical velocities are nonexistent at the surface 
of the earth and the outer boundary of the atmos- 
phere. 

One form of the equation of continuity is 

^ + V.Vp+pV-V+^( p W)=0 (Bl) 



AM/ 10 FT 
56 46 36 31 26 21 16 11 




.0001 



10 5 10 6 

DISTANCE FT 



Figure 15. Nomogram of dj= (aj — aj-i) (KMM- 1 ). 



where p = density of the air, V= horizontal wind 

velocity vector, v= horizontal gradient operator, 

Z= vertical coordinate, W= vertical wind velocity. 

Integration of (Bl) from the LND to the outer 

atmosphere results in 

(pW) LND =f a V'VpdZ+T QdZ+T P V-VdZ. 

J LND jLND<->r ,/LND 

(B2) 

It is usually observed that dp/dt<^py-V so that under 
most circumstances (B2-b) may be disregarded. It 
is also observed that cold air advection below the 
tropopause is accompanied by warm air advection 
in the stratosphere. For such a case (B2-a)— >0. 
With these assumptions (B2) reduces to 



<Wh 



if 

P Jli 



pV-VdZ. 



(B3) 



P J LND 

Horizontal velocity divergence is defined as 

bs drt 

in the natural coordinate system. Evaluation di- 
rectly from the wind field is difficult to accomplish 
with any degree of accuracy. For this reason an 
alternate expression for the integrand of (B3) is 
sought. 
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The vorticity equation is expressed 



r 



dt 



(B4) 



where f is the vertical component of the absolute 
vorticity. The substantial derivative equated to its 
partial components is 



dt dt^ ds T bZ 



(B5) 



Substituting eq 



where s is along the streamlines. 
(B5) in (B3) we have 



(B6) 



The first integral on the right becomes proportion- 
ately smaller with increasing height above the LND 
and at the 300-mb level is generally quite small com- 
pared with the second integral. This is merely a 
statement that, in general, the winds blow througli 
the pressure contours faster than the pressure con- 
tours are moving. The third integral is small because 
in the integral over the entire upper atmosphere the 
only contribution can pass through the lower bound- 
ary. Because the vertical velocity is several orders 
of magnitude smaller than horizontal components, 
the contribution cannot be great. 
We can then write 



(WO, 



PJLND f 



V^-dZ- 
ds 



(B7) 



The absolute vorticity, f=/+f r , and its gradient 
df/bs = 5//ds + dfr/^ may be evaluated directly from 
a constant pressure chart. The gradient of the 
Coriolis term, df/ds is, under most circumstances, at 
least an order of magnitude less than the gradient of 
relative vorticity, df r /ds. 

The magnitude of the vertical component of relative 
vorticity is r r =(dt>/d&) — (bu/dy) Substituting in eq 
(B7) we have finally 



<W)i 



pJlnd f i>s \bx by) 



dZ (B8) 



If the assumptions stated at first are reasonable 
then the sign of 

(W) LND =-sign of [| (Irfj)! >>»W 

Further, the 300-mb level is near the Jetstream 
(region of maximum winds) and (with somewhat 
less certainty) one can also say 

relative magnitude (Wlnd= relative magnitude of 



fl(s-S) 



) 



7.2. Computation of f r 

The computation of relative vorticity on a con- 
stant pressure chart is greatly simplified if one makes 
use of the geostrophic wind approximation 



tt — V dh 

Uo ~~fW 



tt__9 [dh 
r "7 dx 



The finite difference equation for the relative vortic- 
ity is then 



f, 



-ft 



Ah 



(Ax) 2 



Ah \ 

'(Ayr) 



(B9) 



The vorticity field is computed by reading the con- 
stant pressure surface height at evenly spaced grid 
points and carrying out the computation 

where/ is the Coriolis parameter (twice the angular 
velocity of the earth times the sine of the latitude), 
g the acceleration of gravity, D is the grid spacing, 
and Hi, H 2 , H 3 , Hi, H 5 are the pressure surface 
heights at the points P 1 ... P 5 as shown in figure 
16. The computed vorticity is the average value for 
the grid space and is plotted at P 5 . 




Figuke 16. Schematic of grid used to compute relative vorticity 
on a constant pressure surface. 



In this report the relative vorticity was computed 
using a grid size of 5° latitude for a grid centered over 
Catalina Island and another grid 5° latitude upwind. 
The vorticity gradient was assumed to be the finite 
difference between f r at the two positions divided 
by the 5° latitude distance between points. 
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